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Following a brief introduction to early Universe cosmology, the current of status 
of big bang nucleosynthesis is reviewed and the concordance between theory and 
observation is examined in detail. The abundances of ''He and '^Li determine the 
value of the baryon-to-photon ratio, r; to be relatively low, r; 1.8 X 10"^'', and 
agrees with some recent measurements of D/H in quasar absorption systems. These 
results have far reaching consequences for galactic chemical evolution, the amount 
of baryonic dark matter in the Universe and on the allowed number of degrees of 
freedom in the early Universe. Finally, motivations for cosmological dark matter 
will be reviewed with special emphasis placed on supersymmetric candidates. 



1 Introduction 

To set the framework for the subsequent discussion on Big Bang Nucleosyn- 
thesis (BBN) and dark matter, it will be useful to briefly review some general 
aspects of the standard hot big bang model. BBN occurs at the relatively 
late time of order 1 minute after the big bang, however, the initial conditions 
which determine BBN, were set up much earlier. These include, an isotropic 
background, a finite baryon density and a period of equilibrium. The standard 
big bang model assumes homogeneity and isotropy, though it is reasoned that 
these conditions were achieved after a period of inflation. Given homogeneity 
and isotropy, space-time can be described by the Friedmann- Robertson- Walker 
metric which in co-moving coordinates is given by 



ds2 ^ -dt^ +R^{t) 



dr 
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where R{t) is the cosmological scale factor and k is the three-space curvature 
constant (fc = 0, +1,— 1 for a spatially flat, closed or open Universe), k and 
R are the only two quantities in the metric which distinguish it from flat 
Minkowski space. It is also common to assume the perfect fluid form for 
the energy-momentum tensor described by, the isotropic pressure, p, and the 
energy density p. Einstein's equations yield the Friedmann equation, 
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and 

UnGNip + Sp) (3) 



where A is the cosmological constant, or equivalently from the conservation of 
the energy momentum tensor 

P = ~iH{p + p) (4) 

These equations form the basis of the standard big bang model. 

At early times (t < 10^ yrs) the Universe is thought to have been dom- 
inated by radiation so that the equation of state can be given hy p ~ p/3. 
Neglecting the contributions to H from k and A (this is always a good approx- 
imation for small enough R), we find that 

R{t) - 

and p ^ i?^* so that t ^ (3/327rG'Arp)^/^. Similarly for a matter or dust 
dominated Universe with p = 0, 

R{t) - ^2/3 (g) 

and p ^ i?^"^. The Universe makes the transition between radiation and matter 
domination when p^ad = Pmatter or when T ~ few x 10'^ K. 

In the absence of a cosmological constant, one can define a critical energy 
density pc such that p = pc for fc = 

Pc ^ 3H^/8t:Gn (7) 

In terms of the present value of the Hubble parameter this is, 

Pc = 1.88 X IQ-^^ho^gcm-^ (8) 

where 

ho = iJo/(100kmMpc-^s-i) (9) 
The cosmological density parameter is then defined by 

17=— (10) 

Pc 

in terms of which the Friedmann equation, Eq. (^), can be rewritten as 

in-l)H'^^ (11) 



so that fc = 0, +1, —1 corresnonds to il = 1, 17 > 1 and fl < 1. Broad observa- 
tional limits on ho and fl area 



0.4 <ho< 1.0 



0.1 < f7 < 2 



(12) 



The value of fl, at least on relatively small scales, seems to depend on 
scale. Indeed, the contribution to Q from visible matter associated with stars 
and hot gas is quite small, fl « 0.003 — 0.01. On somewhat larger scales, that 
of galactic halos or small groups of galaxies, « 0.02 — 0.1. On galaxy cluster 
scales, it appears that may be asJarge as 0.3. And while there is some 
evidence (see the lectures of N. Kaisen in these proceedings), the observations 
are far from conclusive in indicating a value of as large as 1. It is however 
possible to obtain a bound on the product, flh^ from 



(for A = 0). For to > 12Gyr, and 17 < 1, nh'^ < 0.3 (This is true even if 
A^O.) 



As indicated above, BBN takes place during the radiation dominated epoch 
which lasts roughly to the period of recombination (somewhat earlier when 
dark matter is included) which occurs when electrons and protons form neutral 
hydrogen through + p ^ H +7 at a temperature Tr^ few xlO"^ K ^1 
eV. For T < Tfi, photons are decoupled while for T > Tji, photons are in 
thermal equilibrium. Today, the content of the microwave background consists 
of photons with Tq ~ 2.728 ± .002 Ku. The energy density of photons in the 
background can be calculated from 



and g-y = 2 is the number of spin polarizations for the photon, E-.^ = q is just 
the photon energy (momentum). (I am using units such that % ^ c = = ^ 
and will do so through the remainder of these lectures.) Integrating (|lj) gives 




(13) 




(14) 



where the density of states is given by 



dn~. — -—^[exp{Ej/T) — 1] ^q^dq 



(15) 



(16) 



which is the familiar blackbody result. 
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In general, at very early times, at very high temperatures, other particle 
degrees of freedom join the radiation background when T^rrii for each particle 
type i if that type is brought into thermal equilibrium through interactions. 
In equilibrium the energy density of a particle type i is given by 

Edriq^ (17) 



and 

9i r u ,, ^ /T^i _i_ ii-i„2 

27r2 



driq^ = ^[^xp[{Eq, - ^^^)|T] ± ir^q^dq (18) 



where again gi counts the total number of degrees of freedom for type i, 

/ii is the chemical potential if present and ± corresponds to either Fermi or 
Bose statistics. 

In the limit that T ^ the total energy density can be conveniently 
expressed by 

(20) 

where gB{F) are the total number of boson (fermion) degrees of freedom and 
the sum runs over all boson (fermion) states with to <C T. The factor of 7/8 



is due to the difference between the Fermi and Bose integrals. Equation (20) 
defines N(T) by taking into account new particle degrees of freedom as the 
temperature is raised. 

In the radiation dominated epoch, eq. (^ can be integrated (neglecting the 
T-dependence of TV) giving us a relationship between the age of the Universe 
and its temperature 
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T-^ (21) 



^327r3GjviV(r)^ 

Put into a more convenient form 

tT!,^y - 2.4[iV(r)]-i/2 (22) 

where t is measured in seconds and TmcV in units of MeV. 

The value of N(T) at any given temperature depends on the particle 
physics model. In the standard SU{3) x SU{2) x U{1) model, we can specify 
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Table 1; Effective numbers of degrees of freedom in the standard model. 



Temperature 


New Particles 


AN{T) 


T < me 


7's + i/'s 


29 


rrie < T < 


e± 


43 


<T < rriTT 




57 


m.„ <T < Tc* 


tt's 


69 


Tc ^ T K. ?Tistrange 


- tt's + u,u,d,d + gluons 


205 


TTls <^ T <i TTlcharm 


s, s 


247 


rUc <T < nir 


c, c 


289 


m.r <T < mbottom 


r± 


303 


mi, <T < mw,z 


b,b 


345 


'mw,z <T < mtop 




381 


mt <T < mniggs 


t,t 


423 


Mh <T 


H° 


427 



*Tc corresponds to the confinement-deconfinenient transition between quarks 
and liadrous. N{T) is shown in Figure 1 for Tc = 150 and 400 MeV. It has 
been assumed that muiggs > mtop- 

N(T) up to temperatures of 0(100) GcV. The change in N can be seen in the 
following table. 

At higher tcmpcratTircs (T ^ 100 GcV), N{T) will be model dependent. 
For example, in the minimal SU{5) model, one needs to add to N(T), 24 states 
for the X and Y gauge bosons, another 24 from the adjoint Higgs, and another 
6 (in addition to the 4 already counted in W^, Z and H) from the 5 of Higgs. 
Hence for T > Mx in minimal SU{5), N{T) = 160.75. In a supersymmetric 
model this would at least double, with some changes possibly necessary in the 
table if the lightest supersymmetric particle has a mass below Mh- 

The presence of a particle species in the thermal background assumes ther- 
mal equilibrium and hence interaction rates which are sufiic;iently fast com- 
pared with the expansion rate of the Universe. Roughly, this translates to the 
condition for each particle type i, that some rate F, involving that type be 
larger than the expansion rate of the Universe or 

Ti>H (23) 

in order to be in thermal equilibrium. 

Examples of a processes in equilibrium at early times which drops out of 
equilibrium or decouples at later times are the processes which involve neutri- 
nos. If we consider the standard neutral or charged-current interactions such 
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Log(T/MeV) 

Figure 1: The eflfective numbers of relativistic degrees of freedom as a function of tempera- 
ture. 



ase+ + e ^ly + vore + v-i-^e + u etc., the rates for these processes can be 
approximated by 

r = n{av) (24) 
where (av) is the thermally averaged weak interaction cross section 

{av)^ 0(10-2)TVM,t (25) 

and n is the number density of leptons. Hence the rate for these interactions 
is 

r^k- o(io-2)tVm^ (26) 

The expansion rate, on the other hand, is just 



H = (^^11^^ = (^^N{T)^ ^'\^/Mp ~ lMN{Tf'^T^/Mp. (27) 



1/2 



The Planck mass Mp = G^^'"^ = 1.22 x lO^^ GeV. 
Neutrinos will be in equilibrium when F^k > H on 

T > (500M^)/Mp)i/3~ IMeV. 
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(28) 



The temperature at which these rates are equal is commonly referred to as the 
decoupling or freeze-out temperature and is defined by 

T{Td) = H{Td) (29) 

For temperatures T > T^, neutrinos will be in equilibrium, while for T < 
they will not. Basically, in terms of their interactions, the expansion rate is 
just too fast and they never "see" the rest of the matter in the Universe (or 
themselves). Their momenta will simply redshift and their effective tempera- 
ture (the shape of their momenta distribution is not changed from that of a 
blackbody) will simply fall with T~l/i?. 

Soon after decoupling the pairs in the thermal background begin to 
annihilate (when T < mg). Because the neutrinos arc decoupled, the energy 
released heats up the photon background relative to the neutrinos. The change 
in the photon temperature can be easily computed from entropy conservation. 
The neutrino entropy must be conserved separately from the entropy of inter- 
acting particles. If we denote T>, the temperature of photons, and before 
annihilation, we also have = T> as well. The entropy density of the inter- 
acting particles at T = T> is just 

,^^!^^1(2 + I)(^)r^ (30) 

while at T = T< , the temperature of the photons just after annihilation, 
the entropy density is 

.< = 1^ = 1(2)(^)T| (31) 

and by conservation of entropy .s< = .s> and 

(T</r>)3 = 11/4 (32) 

Thus, the photon background is at higher temperature than the neutrinos 
because the annihilation energy could not be shared among the neutrinos, 
and 

= (4/11)1/% _ 9^ (-33^ 

As we will see below, standard BBN depends on one single parameter, the 
baryon-to-photon ratio, rj = nB/n^- In fact this quantity is really related to 
the net baryon density ub — n§. It would seem however, that today, = 0, 
and that nBjn^ ~ 2 x 10"^°. In the absence of baryon number violation, this 
ratio is roughly constant (in an adiabatically expanding Universe, the ratio 
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of the baryon density to entropy density is constant). The question we face 
therefore, is what determines these initial conditions, rather than say, ns — ng 
and/or ns/n^ = 1- 

Let us for the moment assume that in fact rj — 0. We can compute the 
final number density of nucleons left over after annihilations of baryons and 
antibaryons have frozen out. At very high temperatures, T > 1 GeV (but 
neglecting the quark- hadron transition), nucleons were in thermal equilibrium 
with the photon background and ns ^ ng ~ (3/2)n-y (a factor of 2 accounts 
for neutrons and protons and the factor 3/4 for the difference between fermi 
and bose statistics). As the temperature fell below niN, annihilations kept the 
nuclcon density at its equilibrium value iriB/n^) — {niN /T)^^'^exp{—mN /T) 
until the annihilation rate Ta — nB'm~^ fell below the expansion rate. This 
occurred at T ~ 20 McV. However, at this time the nuclcon number density 
had already dropped to 

ns/n^ = 71^/% ~ 10"^^ (34) 

which is eight orders of magnitude too small in addition to the problem of 
having to separate the baryons from the antibaryons. If any separation did 
occur at higher temperatures (so that annihilations were as yet incomplete) 
the maximum distance scale on which separation could occur is the causal 
scale related to the age of the Universe at that time. At T = 20 MeV, the 
age of the Universe was only t = 2 x 10~^ sec. At that time, a causal region 
(with distance scale defined by 2ct) could only have contained 10~^Mq which 
is very far from the galactic mass scales on which we are asking for separations 
to occur, IO^^Mq. In spite of all of these problems, 77 = implies that the 
Universe as a whole is baryon symmetric, thus unless baryons are separated 
on extremely large (inflationary) domains, in which case we might just as well 
worry again about ry 7^ 0, there should be antimatter elsewhere in the Universe. 
To date, the only antimatter observed is the result of a high energy collision, 
either in an accelerator or in a cosmic-ray collision in the atmosphere. There 
has been no sign to date of any primary antimatter, such as an anti-helium 
nucleus a found in cosmic-rays. 

The production of a net baryon asymmetry requires baryon number violat- 
ing interactions, C and CP violation and a departure from thermal equilibriunxl. 
The first two of these ingredients are contained in Grand Unified Theories 
(GUTs), the third can be realized in an expanding universe where as we have 
seen, it is not uncommon that interactions come in and out of equilibrium. In 
SU(5), the fact that quarks and leptons are in the same multiplets allows for 
baryon non-conserving interactions such as e~ + d u + u, etc., or decays of 
the supermassive gauge bosons X and Y such as X ^ + d,u + u. Although 
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today these interactions are very ineffective because of the very large masses 
of the X and Y bosons, in the early Universe when T ^ A4x ^ 10^^ GcV these 
types of interactions should have been very important. C and CP violation is 
very model dependent. In the minimal SU(5) model the magnitude of C and 
CP violation is too small to yield a useful value of rj. The C and CP violation 
in general comes from the interference between tree level decay diagrams and 
their one loop corrections. 

The departure from equilibrium is very common in the early Universe when 
interaction rates cannot keep up with the expansion rate. In fact, the simplest 
(and most useful) scenario for baryon production makes use of the fact that 
a single decay rate goes out of equilibrium. It is commonly referred to as the 
out of equilibrium decay scenario eI. The basic idea is that the gauge bosons X 
and Y (or Higgs bosons) may have a lifetime long enough to insure that the 
inverse decays have already ceased so that the baryon number is produced by 
their free decays. 

More specifically, let us call X either the gauge boson or Higgs boson, 
which produces the baryon asymmetry through decays. Let a be its coupling 
to fermions. For X a gauge boson, a will be the GUT fine structure constant, 
while for X a Higgs boson, (47ra)^/^ is its Yukawa coupling to fermions. The 
decay rate for X will be 

Yd ^ aMx (35) 

However decays can only begin occurring when the age of the Universe is longer 
than the X lifetime P^^, i.e., when P^ > H 

aMx > N{Tf''^T^lMp (36) 

or at a temperature 

<, aMxMpN{T)-^/^. (37) 

Scatterings on the other hand proceed at a rate Pg ~ a^T^/Mj^ and hence are 
not effective at lower temperatures. To be in equilibrium, decays must have 
been effective as T fell below Mx in order to track the equilibrium density 
of X's (and X's). Therefore, the out-of-equilibrium condition is that at T = 
Mx,Td <H or 

Mx > aMp{N{Mx)y^''^ ^ lO^^aCeV (38) 

In this case, we would expect a maximal net baryon asymmetry to be produced. 

In the out-of-equilibrium decay scenario □, the total baryon asymmetry 
produced per X, X pair is proportional to AB = (f — r) as can easily be seen 
from Figure 2. Here r (f) is the relevant branching ratio for the decay of X 
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(X). C and CP violation translate into the condition that r ^ f. If decays 
occur out-of-equilibrium, then at the time of decay, nx ~ ai T < Mx- We 
then have 

Some comments regarding the baryogenesis mechanism described above 
are in order: Initial conditions are not important. Equilibrium prior to decay 
erases any net baryon number. The out of equilibrium decays of X and Y then 
generate a new asymmetry. Therefore, the final resulting asymmetry is always 
the same, given by say, (|9|). While ns/s ~ 10"^" is possible, the final result 
is strongly model dependent. Lshould add that several other mechanisms for 
baryogenesis are known to existQ, but all require the same three key ingredients. 

There are several long standing cosmological problems whose solution we 
associate with inflation (see refs. 2 and 8). Here I will describe only one of 
these. The curvature problem (or flatness problem or age problem) can mani- 
fest itself in several ways. For a radiation dominated gas, the entropy density 
s ~ T'^ and conservation of entropy implies R ~ T"^. Thus assuming an 
adiabatically expanding Universe, the quantity k = k/ R^T^ is a dimensionless 



constant. If we now apply the limit in Eq. (12) to Eq. (hlf) we flnd 



(Qo - l)Ho 



2 



R^T^ To 
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< 2 X 10"^*^ (40) 



This limit on k represents an initial condition on the cosmological model. The 
problem is: what physical processes in the early Universe produced a value of 
k so extraordinarily close to zero (or fl close to one)? A more natural initial 
condition might have been fc ~ 0(1). In this case the Universe would have 
become curvature dominated at T ~ 10~"'^Afp. 

It is important to note that il is a function of time or of the scale factor. 
The evolution of is shown in Figure 3 for A = 0. For a spatially flat Universe, 
f2 = 1 always. When k = +1, there is a maximum value for the scale factor R. 
At early times (small values of i?), fl always tends to one. Note that the fact 
that we do not yet know the sign of fc, or equivalently whether f2 is larger than 
or smaller than unity, implies that we are at present still at the very left in the 
figure. What makes this peculiar is that one would normally expect that the 
sign of k to become apparent after a Planck time of lO^''^ s. It is extremely 
puzzling that some 10^° Planck times later, we still do not know the sign of fc. 

The curvature problem discussed above among others can be. neatly re- 
solved if the Universe underwent a period of cosmological inflationQ. During a 
phase transition, our assumptions of an adiabatically expanding universe may 
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Figure 3: The evolution of the cosmological density parameter, Q, as a function of the scale 
factor for a closed, open and spatially flat Universe. 

not be valid. If we look at a scalar potential describing a phase transition from 
a symmetric false vacuum state (ip) — for some scalar field (p to the broken 
true vacuum at (</>) = v, and suppose we find that upon solving the equations 
of motion for the scalar field that the field evolves slowly from the symmetric 
state to the global minimum (this will depend on the details of the potential). 
If the evolution is slow enough, the universe may become dominated by the 
vacuum energy density associated with the potential near (/> sa 0. The energy 
density of the symmetric vacuum, V(0) acts as a cosmological constant with 

A = 8nV{0)Mp^ (41) 

During this period of slow evolution, the energy density due, to say, radiation 
will fall below the vacuum energy density, p ^ V{Q). When this happens, the 
expansion rate will be dominated by the constant V(0) and from Eq. (^ we 
find an exponentially expanding solution 

R{t) - * (42) 

When the field evolves towards the global minimum it will begin to oscillate 
about the minimum, energy will be released during its decay and a hot thermal 
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universe will be restored. If released fast enough, it will produce radiation at a 
temperature Tr^ < V{0). In this reheating process entropy has been created 
and {RT)f > {RT)i. Thus we see that during a phase transition the relation 
RT ^ constant, need not hold true and thus our dimensionless constant k may 
actually not have been constant. 

If during the phase transition, the value of RT changed by a factor of 
0(10^^), the cosmological problems would be solved. The isotropy would in a 
sense be generated by the immense expansion; one small causal region could 
get blown up and hence our entire visible Universe would have been at one 
time in thermal contact. In addition, the parameter k could have started out 
0(1) and have been driven small by the expansion. 

2 Big Bang Nucleosynthesis 

The concordance between big bang nucleosynthesis (BBN) theory and obser- 
vation has been the subject of considerable recent debate. It is clear however, 
that the real questions lie not with the concordance between BBN and the 
observational data, but rather between the theories of chemical and stellar 
evolution and the data. BBN theory (see for example ref. 9) is quite stable 
in the sense that over time very little in the fundamental theory has changed. 
Cross-sections are now more accurately measured, the neutron mean life is 
known with a much higher degree of precision, and if we restrict our attention 
to the standard model, the number of neutrinos has also been determined. In 
contrast, the status of the observational data has changed significantly in the 
last several years. There is better data on ''He, more data on ^Li, and data 
on D and ^He that was simply non-existent several years ago. For the most 
part, the inferred abundances of "^He and ''Li have remained relatively fixed, 
giving us a higher degree of confidence in the assumed primordial abundances 
of these isotopes as is reflected in their observational uncertainties. Indeed, the 
abundances of '*He and ^Li alone are sufRcient in order to probe and test the 
eory and determine the single remaining parameter in the standard model 
, namely, the baryon-to-photon ratio, rj. In contrast, D and ■^He are highly 
dependent on models of chemical evolution ('^He is in addition dependent on 
the uncertain stellar yields of this isotope). New data from quasar absorp- 
tion systems, on what may be primordial D/H is at this time disconcordant, 
different measurements give different abundances. As a consequence of the 
uncertainties in D and ^He, one can use the predictions based on ^He and ^Li 
in order to construct models of galactic chemical evolution. These results also 
have important implications for the amount of (non)-baryonic dark matter in 
the galaxy and on the number of allowed relativistic degrees of freedom at the 



13 



time of BBN, commonly parameterized as N^. The basic results of BBN are 
summarized in Figure 4, showing the range in 77 where consistency with the 
observations is achieved. 

There has always been an intimate connection between BBN and the mi- 
crowave background as key test to the standard big bang model. Indeed, it 
was the formulation q£ BBN which predicted the existence of the microwave 
background radiationliil. The argument is rather simple. BBN requires temper- 
atures greater than 100 keV, which according to eqs. ( pT| ) and (|2|) corresponds 
to timescales less than about 200 s. The typical cross section for the first link 
in the nucleosynthetic chain is 

(Tv{p + n^D + -i)~bx 10~^"cmVs (43) 

This implies that it was necessary to achieve a density 

nr^ — lO^^cm"^ (44) 

Now the density in baryons today is known approximately from the density of 
visible matter to be ubo ^ 10^^ cm^'^ and since we know that that the density 
n scales as R^^ ^ T^, the temperature today must be 

To = (nBo/")'/'TBBN - lOK (45) 

A pretty good estimate. 

Despite its simplicity, BBN was criticized early on, due to its shortcom- 
ings in being able to produce the observed abundances of all of the element 
isotopes. Attention was therefore turned to stellar nucleosynthesis. However, 
while the elements from helium on up can be and are produced in stars, no 
other astrophysical site has ever survived for the production of deuterium. In 
addition, if one assumes that "'He is entirely of stellar origin, one should be able 
to find places in the Universe in which the "^He mass fraction is substantially 
below 25%. The ^He data shown in Figure 5, emphasizes the fact that indeed 
no such region with low ^He has ever been observed and that (together with 
the need to produce D) leads one to conclude that BBN nucleosynthesis is a 
necessary component in any cosmological model. 

Before commencing with the direct comparison between theory and ob- 
servations, it will be useful to briefly review the main events leading to the 
synthesis of the light elements. Conditions for the synthesis of the light ele- 
ments were attained in the early Universe at temperatures T < 1 MeV. At 
somewhat higher temperatures, weak interaction rates were in equilibrium, 
thus fixing the ratio of number densities of neutrons to protons. At T ^ 1 
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Figure 5: The helium (Y) vs oxygen (O/H) abundances in extragalactic HII regions em- 
phasizing the lack of low *He regions. 

MeV, {n/p) ~ 1. As the temperature fell and approached the point where 
the weak interaction rates were no longer fast enough to maintain equilibrium, 
the neutron to proton ratio was given approximately by the Boltzmann factor, 
(n/p) ~ e"'^™/^, where Am is the neutron-proton mass difference. The final 
abundance of "^He is very sensitive to the (n/p) ratio. 

The nucleosynthesis chain begins with the formation of deuterium through 
the process, p + n —^ Y) -I-7. However, because the large number of photons rel- 
ative to nucleons, ri~^ — n^/nB ~ 10^°, deuterium production is delayed past 
the point where the temperature has fallen below the deuterium binding en- 
ergy, Eb = 2.2 MeV (the average photon energy in a blackbody is ~ 2.7T). 
When the quantity 77~^exp(— £^b/T) ^ 1 the rate for deuterium destruction 
(D -I-7 — > p -t- n) finally falls below the deuterium production rate and the 
nuclear chain begins at a tempiecature T ~ O.lMeV. The nuclear chain in 
BBN calculations was extended 113 and is shown in Figure 6. 

The dominant product of big bang nucleosynthesis is "^He resulting in an 
abundance of close to 25% by mass. This quantity is easily estimated by 
counting the number of neutrons present when nucleosynthesis begins. When 
the weak interaction rates responsible for n — p equilibrium freeze-out, at T « 
0.8 MeV, the neutron to proton ratio is about 1/6. When free neutron decays 
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Figure 6: The nuclear network used in BBN calculations. 
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prior to deuterium formation are taken into account, the ratio drops to (n/p) ~ 
1/7. Then simple counting yields a *He mass fraction 



^ « 0.25 (46) 

[1 + (n/p)] 

In the standard model, the ^He mass fraction depends primarily on the 
baryon to photon ratio, 77 as it is this quantity which determines the onset of 
nucleosynthesis via deuterium production. For larger values of 77, the quantity 
ri~^exp{—EB/T) ^ 1 is smaller, and hence the nuclear chain may begin at 
a higher temperature. As a result, the (n/p) ratio is higher, producing more 
^He. But because the (n/p) ratio is only weakly dependent on 77, the ^He mass 
fraction is relatively flat as a function of 77. When we go beyond the standard 
model, the ^He abundance is very sensitive to changes in the expansion rate 
which can be related to the effective number of neutrino flavors as will be 
discussed below. Lesser amounts of the other light elements are produced: D 
and '^He at the level of about 10~^ by number, and ^Li at the level of_Lp~^° 
by number. These abundances (along with ^Li) are shown in Figure 7ll3. In 
Figure 8, the produced abundances of the intermediate mass isotopes ^Be, ^°B, 
^^B are also shown. These abundances are far below the observed values and 
it is believed that these isotopes are formed in cosmic ray nucleosynthesis. 

For the comparison with the observations, I will use the resulting abun- 
dances of the light elements shown in Figure 4. The curves for the *He mass 
fraction, Y, bracket the computed range_based on the uncertainty of the neu- 
tron mean-life which has been taken asllj r„ = 887 ± 2 s. Uncertainties in the 
produced ''Li abundances have been adopted from the results in Hata et al. 
Uncertainties in D and '^He production are small on the scale of this figure. 
The boxes correspond to the observed abundances and will be discussed below. 



2.1 Light Element Abundances 

There is now a good collection of abundance information on the ^He mass frac- 
tinn V" 0/H, and N/H in over 70 extragalactic HII (ionized hydrogen) regions 
E3il3'tZl. The ^He abundance in very low metallicity regions is best determined 
from observations of Hell — > Hel recombination lines. The observation of the 
heavy elements is important as the helium mass fraction observed in these HII 
regions has been augmented by some stellar processing, the degree to which 
is given by the pxjiigen and nitrogen_abundances. In an extensive study based 
on the data in Ej^EEl, it was found tB that the data is well represented by a 
linear correlation for Y vs. 0/H and Y vs. N/H. It is then expected that 
the primordial abundance of *He can be determined from the intercept of that 
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Figure 7: The light element abundances from big bang nucleosynthesis as a function of rj, 
including ®Li. 
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Figure 8: The intermediate mass element abundances from big bang nucleosynthesis as a 
function of r]. 
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relation. The overall result of that analysis indicated a primordial mass frac- 
tion, Yp — 0.232 ± jli)03. The stability of this fit was verified by a statistical 
bootstrap analysis Ej showing that the fits were not overly sensitive to any 
particular HII region. 

To make use of the ''^He data, it is crucial to obtain high quality and 
very low metallicity data. In principle, any one HII region (with non-zero 
metallicity) should provide an upper limit to Yp since some stellar processing 
has taken place augmenting the primordial value. Thus the determination of 
Yp by an extrapolation to zero metallicity could be avoided by the observations 
of either low metallicity or low helium HII regions. For a very low metallicity 
HII region such an upper limit may even provide a reasonable estimate of Yp. 
For example, I Zw 18, is the lowest metallicity extragalactic HII region with 
''He observed. In I Zw 18, there are 2 separate regions, and there have been 
five distinct measurements of ''He. In the HII regions designated as NW and 
SE, the five measured abundances are shown in the table below. This data 



Table 1: Data for I Zw 18 
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.012 


ref. 
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.231 ± 
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ref. 
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NW 


.232 ± 


.008 


ref. 
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SE 


.231 ± 


.012 


ref. 


15 


SE 


.230 ± 


.008 


ref. 
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leads to a total average for I Zw 18 of YJ, = 0.230 ± 0.004. There are several 
other regions with low *He abundances and thepjalue of Y derived from the 
average of several HII regions is also of interest If we successively include 
additional HII regions with higher Y, as more regions are included, the mean 
value (weighted) of Y will increase, but if the errors are statistical, the error 
in the mean will decrease. As a result, for N HII regions the one-(or two-)a 
upper bound to Y will first decrease with N, then level off and, as N is further 
increased, it will eventually increase monotonically. This behavior is seen in 
Figure 9 where weighted means are shown, and the Icr bounds to the weighted 
means of Y derived from the N lowest helium abundance HII regions. Note 
that for 2 < iV < 13, the mean varies from 0.229 to 0.231 while for 2 < iV < 14, 
(y) < 0.236(2ct). It is not unreasonable to infer from these results that, 

Yp < 0.230 ± 0.003 (47) 

with, Yp"' < 0.236. If, instead, we take the weighted means of the regions 
with the lowest values of 0/H, we obtain a similar result. This illustrates 
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Figure 9: The running average (weighted means) of the *He abundance, Y, for the first N 
(lowest Y) HIT regions (from ref. 20). Also shown arc the la bounds to the weighted means. 



the potentially great value of very careful analyses of a handful of the lowest 
metallicity (lowest Y) HII regions. 

Although the above estimates on Yp are consistent with those based on a 
linear extrapolation of the data in refs. 15 and 16, it has been claimed, that 
the new data in ref. 17, leads to a significantly higher value for Yp (in excess 
of 24%). However, in a recent analysis EJ it has been shown that the data in 
ref. 17 is entirely consistent with the data in refs. 15 and 16 as can be seen in 
Figure 10 where the Y versus 0/H data used in ref. 18 (open circles, based on 
the data in refs. 15 and 16)) is shown along with the new data (filled circles, 
from ref. 17). The crossed circles are the 10 HII regions that Izotov et al. 
excluded from their own analyses. Note that where there is overlap in 0/H, 
the newer Y values are intermingled with those from refs 15 and 16. It was 
arguedc3 that because the data in ref. 17 does not extend down to sufficiently 
low metallicity, taken alone it can not be used to determine Yp. It was also 
shownE3 that one can not argue for a systematic shift in Y on the basis of new 
atomic data calculations as was done in ref. 17. On the other hand, this data 
may be jCombined with the previous data in which case one finds a "^He mass 
fractioncj based on 62 distinct HII regions 

Yp = 0.234 ± 0.002 ± 0.005 (48) 
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Figure 10: The helium (Y) and oxygen (O/H) abundances from ref. 20 in extragalactic 
HII regions, (open circles) from refs. 15 and 16, and from ref. 17 (filled circles). Regions 
excluded by ref. 17 are shown as crossed circles. Lines connect the same regions observed 
by different groups. 
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The first uncertainty is purely statistical and the second uncertainty is an esti- 
mate of the systematic uncertainty in the primordial abundance determination. 
The solid box for ''He in Figure 4 represents the range (at 2(Tstat) from (^8|). 
The dashed box extends this by Csys ■ A somewhat lower primordial abundance 
Yp — 0.230 ± .003 is found by restricting to the 32 most metal poor regions 

The ^Li abundance is also reasonably well known. In old, hot, population- 

11 stars, ^Li is found to have a very nearly uniform abundance. For stars with 
a surface temperature T > 5500 K and a metallicity less than about l/20th 
solar (so that effects such as stellar convection may not be important), the 
abundances show little or no dispersion beyond that which is consistent with 
the errors of individual measurements. Indeed, as detailed in ref. 21, much of 
the work concerning ^Li has to do with the presence or absence of dispersion 
and whether or not there is in fact some tiny slope to a [Li] = log ^Li/H + 

12 vs. T or [Li] vs. [Fe/H] relationship ([Fe/H] is the log of the Fe/H ratio 
relative to the solar value). 

There is ^Li data from nearly 100 halo stars, from a variety of sources. 
When the Li data from stars with [Fe/H] < -1.3 is plotted as a function of 
surface temperature, one sees a plateau emerging for T > 5500 K as shown 
in Figure 11 for the data taken from ref. 22. As one can see from the figure, 
at high temperatures, where the convection zone does not go deep below the 
surface, the Li abundance is uniform. At lower temperatures, the surface 
abundance of Li is depleted as Li passes through the interior of the star and 
is destroyed. The lack of dispersion in the plateau region is evidence that this 
abundance is indeed primordial (or at least very close to it). Another way to 
see the plateau is to plot the Li abundance data as a function of metallicity, this 
time with the restriction that T > 5500 K as seen in Figure 12. In this case, 
the plateau emerges at low [Fe/H] as would be expected. At higher [Fe/H], the 
convection zone remains below the surface only for much hotter stars. Thus, 
for [Fe/H] > -1.3, the effects of depletion are seen. Also apparent in this figure 
is that at higher metallicities there is evidence for the production of Li which 
rises by over an order of magnitude at solar metallicity. 

I will use the value given in ref. 22 as the best estimate for the mean ^Li 
abundance and its statistical uncertainty in halo stars 

Li/H= (1.6 ±0.li°;™)x 10-1" (49) 

The first error is statistical, and the second is a systematic uncertainty that 
covers the range of abundances derived by various methods. The Li abun- 
dance is somewhat sensitive to stellar parameters such as the assumed surface 
temperature, the metallicity and the surface gravity. The greatest model de- 
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Figure 11: The Li abundance in halo stars with [Fe/H] < -1.3, as a function of surface 
temperature 
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Figure 12: The Li abundance in halo stars with T > 5500 K, as a function of mctallicity. 
Filled diamonds represent upper limits. 



25 



pendence is on the conversion of the observed colors (B-V) to temperature. 
For example, in a sample of 55 stars takeii from the papers of ref. 23, one 
finds [Li] = 2.08 ± 0.02. From Thorburn'£J sample of 74 stars one finds [Li] 
= 2.26 ± 0.01. And ref. 22, which comprises 41 stars gives [Li] = 2.21 ± 0.01. 
The Li abundance used in (^) is derived from this last data set. The solid 
box for ''Li in Figure 4 represents the 2(Tstat + osys range from . 

The third set of errors deals with the possible depletion and/or production 
of ^Li. These uncertainties are however limited. As was noted above, the lack of 
dispersion in the Li data limits the amount of depletion. In addition, standard 
stellar modelO predict that any depletion of ^Li would be accompanied by a 
very severe depletion of ^Li. Ifatil recently, ^Li had never been observed in hot 
pop II stars. The observatiorc3 of ^Li (which turns out to be consistent with 
its origin in cosmic-ray nucleosynthesis and with a small amount of depletion 
as expected from standard stellar models) is another good indication that ^Li 
has not been destroyed in these starcZl. 

Aside from the big bang, Li is produced together with Be and B in cosmic 
ray spallation of C,N,0 by protons and a-particles. Li is also produced by 
a — a fusion. Be and B have been observed in these same pop II stars and in 
particular there are a dozen or so stars in which both Be and ^Li have been 
observed. Thus Be (and B though thereJs still a paucity of data) can be used 
as a consistency check on primordial Lit3. Based on the Be abundance found 
in these stars, one can conclude that no more than 10-20% of the ''Li is due to 
cosmic ray nucleosynthesis leaving the remainder (an abundance near 10^ ^*') 
as primordial. The third set of errors in Eq. (^ ) accounts for the possibility 
that as much as half of the primordial ''Li has been destroyed in stars, and that 
as much as 20% of the observed ''Li may have been produced in cosmic ray 
collisions rather than in the Big Bang. The dashed box in Figure 4, accounts 
for this additional uncertainty. For ^Li, the uncertainties are clearly dominated 
by systematic effects. 

Turning to D/H, we have three basic types of abundance information: 1) 
ISM data, 2) solar system information, and perhaps 3) a primordial abundance 
from quasar absorption systems. The best measurement for ISM D/H iscB 

(D/H)isM - 1.60 ± Omt'^ fa x 10"^ (50) 

Because there are no known astrophysical sites for the production of deuterium, 
all observed D must be primordial. As a result, a firm lower bound from 
deuterium establishes an upper bound on i] which is robust and is shown by the 
lower right of the solid box in Figure 4. The solar abundance of D/H is inferred 
from two distinct measurements of '^He. The solar wind measurements of '^He 
as well as the low temperature components of step- wise heating measurements 
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of '^He in meteorites yield the presolar (D + "^Hc) /H ratio, as D was efhciently 
burned to '^He.ixi the Sun's pre-main-sequcncc phase. These measurements 
indicate thatBEl 

= (4.1 ±0.6 ±1.4) X 10"^ (51) 



The high temperatum components in meteorites are believed to yield the true 
solar -^He/H ratio ofE3'EJ 

= (1.5 ±0.2 ±0.3) X 10"^ (52) 



The difference between these two abundances reveals the presolar D/H ratio, 
giving, 

(D/H)0 « (2.6±0.6± 1.4) X 10^^ (53) 

It should be noted that recent measurements of surface abundances Q^flD 
on Jupiter show a somewhat higher value for D/H, D/H = 5 ± 2 x 10^^ e3. If 
this value is confirmed and if fractionation does not significantly alter the D /H 
ratio (as it was suspected to for previous measurements involving CH3D), it 
may have an important impact on galactic chemical evolution models. This 
value is marginally consistent with the inferred meteoritic values. 

Finally, there have been several recent reported measurements of D /H is 
high redshift quasar absorption systems. Such measurements are in principle 
capable of determining the primordial value for D/H and hence rj, because 
of the strong and monotonic dependence of D/H on 77. However, at present, 
detections of D /H using quasar absorption systems do not yield a conclusive 
value for D/H. As such, it should be cautioned that these values may not turn 
out to represent the true primordial value pad it is very unlikely that both are 
primordial and indicate an inhomogeneitytS. The first of these measurements 
O indicated a rather high D/H ratio, D/H « 1.9 - 2.5 xlO~*. Other high 
D/H ratios were reported inH. However, there are reported low values of D/H 
in other such systems with values D/H ~ 2.5 x 10^^, significantly lower 
than the ones quoted above. Though this primordial D/H value is consistent 
with the solar and present values of D/H, it is not consistent (at the 2 a level) 
with the determinations of D/H in Jupiter, if they are correct. It is generally 
regarded that the high D /H measurements are of lesser quality while a recent 
observationLJ of the hydrogen column density in one of low D /H clouds found 
a significantly lower value and claimed a lower limit on D/H of 4 x 10~^. 

The range of quasar absorber D/H is shown by the dashed box in Figure 
4. It is probably premature to use either of these values as the primordial 
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D/H abundance in an analysis of big bang nucleosynthesis, but it is certainly 
encouraging that future observations may soon yield a firm value for D/H. It 
is however important to note that there does seem to be a trend that over the 
history of the Galaxy, the D /H ratio is decreasing, something we expect from 
galactic chemical evolution. Of course the total amount of deuterium astration 
that has occurred is still uncertain, and model dependent. 

There are also several types of '^He measurements. As noted above, me- 
teoritic extractions yield a presolar value for ^He/H as given in Eq. (|5^). Ii 
addition, there are several ISM measurements of '^He in galactic HII regions !f 
which show a wide dispersion which may be indicative of pollution or a biasE; 

~ 1 - 5 X 10"^ (54) 

HII 

There is also a recent ISM measurement of ^^1100 with 

= 2.lti X 10"^ (55) 

ISM 

Finally there are observations of "^He in planetary nebulaeEil which show a very 
high ^He abundance of^He/H ~ 10~^. 

Each of the light element isotopes can be made consistent with theory for 
a specific range in 77. Overall consistency of course requires that the range in 77 
agree among all four light elements. However, as will be argued below D and 
■^He are far more sensitive to chemical evolution than ^He or ^Li and as such 
the direct comparison between the theoretical predictions of the primordial 
abundances of D and '^He with the observational determination of their abun- 
dances is far more difficult. Therefore in what follows I will for the most part 
restrict the comparison between theory and observation to the two isotopes 
who suffer the least from the effects of chemical evolution. 





2.2 Chemical Evolution 

Because we can not directly measure the primordial abundances of any of the 
light element isotopes, we are required to make some assumptions concerning 
the evolution of these isotopes. As has been discussed above, "^He is produced 
in stars along with oxygen and nitrogen. ^Li can be destroyed in stars and pro- 
duced in several (though still uncertain) environments. D is totally destroyed 
in the star formation process and '^He is both produced and destroyed in stars 
with fairly uncertain yields. It is therefore preferable, if possible to observe the 
light element isotopes in a low metallicity environment. Such is the case with 
^He and ^Li, and we can be fairly assured that the abundance determinations 
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of these isotopes are close to primordial. If the quasar absorption system mea- 
surements of D/H stabihze, then this too may be very close to a primordial 
measurement. Otherwise, to match the solar and present abundances of D and 
^He to their primordial values requires a model of galactic chemical evolution. 

The main inputs to chemical evolution models are: 1) The initial mass 
function, 4'{m), indicating the distribution of stellar masses. Typically, a sim- 
ple power law form for the IMF is chosen, (j){m) ^ m^^, with x ~ —2.7. This is 
a fairly good representation of the observed distribution, particularly at larger 
masses. 2) the star formation rate, ijj. Typical choices for a SFR are ip{t) (x a 
or or even a straight exponential e^*/^. a is just the fraction of mass in gas, 
Afgas/Mtot- 3) the presence of infalling or outflowing gas; and of course 4) the 
stellar yields. It is the latter, particularly in the case of '^He, that is the cause 
for so much uncertainty. Chemical evolution models simply set up a series of 
evolution equations which trace desired quantities. For example, a and hence 
the SFR evolve through a relation such as 

^^ = -VW + e(t)+i(t)-o(i) (56) 

where e represents the amount of gas ejected from stars, i is the gas infall rate, 
and o is the gas outflow rate. Similar equations can be developed which trace 
the abundances of the element isotopes c3. 

Deuterium is always a monotonically decreasing function of time in chem- 
ical evolution models. The degree to which D is destroyed, is however a model 
dependent question which depends sensitively on the IMF and SFR. The evo- 
lution of ■^He is however considerably more complicated. Stellar models predict 
that substantial amounts of ■^Hcj-are produced in stars between 1 and 3 M©. 
For M < 8M0, Iben and TruranN calculate 

fHe/H)/ = 1.8 X 10-4 0.7 [(D+ ^HeVH]^ (57) 

so that for example, when t^iq — 3, ((D -f ■^IIe)/H)i = 9 x lO^'^, and the ratio 
of the final abundance of ^He/H to the initial (D + ■^He)/H abundance denoted 
by <?3 is 53(1Mq) = 2.7. The ^He abundance is nearly tripled. It should be 
emphasized that this prediction iiin fact consistent with the observation of 
high ■^He/H in planetary nebulaellJ. 

Generally, implementation of the ■^He yield in Eq. (|57| ) in chemical evolu- 
models leads to an overproduction of ■^He/H particularly at the solar epoch 
This problem is coninounded in models with an intense period of D de- 
struction. In Scully et al. c3, a dynamically generated supernovae wind model 
was coupled to models of galactic chemical evolution with the aim of reducing 
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a primordial D/H abundance of 2 xlO~^ to the present ISM value without 
overproducing heavy elements and remaining within the other observational 
constraints typically imposed on such models. In Figure 13, the evolution of 
D/H and ^He/H is shown as a function of time in several representative models 
with significant deuterium destruction factors (see ref 45. for details). How- 
ever, as one can plainly see, '^He is grossly overproduced (the deuterium data 
is represented by squares and '^He by circles). 

The overproduction of ^He relative to the solar meteoritic value seems to 
be a generic feature of chemical evolution models when ^He production in low 
mass stars is included. This result appears to be independent of the chemical 
evolution model and is directly related to the assumed stellar yields of "^He. 
It has recently been suggested that at least some low mass stars may indeed 
be net destroyers of '^He if one includeaJihe effects of rotational mixing in 
low mass stars on the red giant branch e30. The extra mixing does not take 
place for stars which do not undergo a helium core flash (i.e. stars > 1.7-2 
Mq ). Thus stars with masses less than 1.7 Mq are responsible for the '^He 
destruction. Using the yields of Boothroyd and MalaneyS, it was shown @ 
that these reduced ^He yields in low mass stars can account for the relatively 
low solar and present day ^He/H abundances observed. In fact, in some cases. 
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Figure 14: The evolution of D/H and ^He/H with time using the reduced '^He yields of ref. 
47. The dashed curve is the same as in Figure 13, using standard ^He yields. 



^He was underproduced. To account for the '^He evolution and the fact that 
some low mass stars must be producers of '^He as indicated by the planetary 
nebulae data, it was suggested th. 
(albeit large) of low mass stars 
and "^He/H is shown in Figure 14 



the new yields apply only to a fraction 
The corresponding evolution @ of D/H 



2.3 Likelihood Analysis 

Monte Carlo tcchniqucshaye proven to be a useful form of analysis_for big 
bang nucleosynthesis Ej'Eilo. An analysis of this sort was performedllj using 
only ^He and ^Li. It should be noted that in principle, two elements should be 
sufficient for not only constraining the one parameter (rf) theory of BBN, but 
also for testing for consistency. The procedure begins by establishing likelihood 
functions for the theory and observations. For example, for '^He, the theoretical 
likelihood function takes the form 



L 



BBN 



(r,rBBN)-e-(^"^B-(''))'/2-? 



(58) 



where yBBN(^) is the central value for the "^He mass fraction produced in the big 
bang as predicted by the theory at a given value of t], and ui is the uncertainty 
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Figure 15: Likelihood distribution for each of ^He and '^hi, shown as a function of rj. The 
one-peak structure of the ^He curve corresponds to its monotonic increase with rj, while the 
two peaks for ^Li arise from its passing through a minimum. 

in that value derived from the Monte Carlo calculations and is a measure 
of the theoretical uncertainty in the big bang calculation. Similarly one can 
write down an expression for the observational likelihood function. Assuming 
Gaussian errors, the likelihood function for the observations would take a form 
similar to that in (p8|). 

A total likelihood function for each value of r] is derived by convolving the 
theoretical and observational distributions, which for ^He is given by 



total I, 



(77) = / dYLBBN (Y, Fbbn iri)) Lo{Y, Yo) 



(59) 



An analogous calculation is performed0 for ^Li. The resulting likelihood func- 
tions from the observed abundances given in Eqs. ( p8[ ) and ( ^ ) is shown in 
Figure |l5|. As one can see there is very good agreement between ^He and ''Li in 
the vicinity of 7710 ~ 1.8. The double peaked nature of the ^Li likelihood func- 
tion is due to the presence of a minimum in the predicted lithium abundance. 
For a given observed value of ^Li, there are two likely values of rj. 

The combined likelihood, for fitting both elements simultaneously, is given 
by the product of the two functions in Figure Em and is shown in Figure Eq. 
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Figure 16: Combined likelihood for simultaneously fitting *He and ^Li, as a function of rj. 



From Figure it is clear that ''He overlaps the lower (in 77) ^Li peak, and 
so one expects that there will be concordance in an allowed range of 77 given 
by the overlap region. This is what one finds in Figure [l^, which does show 
concordance and gives a preferred value for 77, ryio = 1.8^^2^ corresponding to 
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(60) 



Thus, we can conclude that the abundances of '*He and ''Li are consistent, 
and select an rjio range which overlaps with (at the 95% CL) the longstanding 
favorite range around 7710 = 3. Furthermore, by finding concordance using 
only ^He and ''Li, we deduce that if there is problem with BBN, it must 
arise from D and ^He and is thus tied to chemical evolution or the stellar 
evolution of "^He. The most model-independent conclusion is that standard 
BBN with Ni, = 3 is not in jeopardy. It is interesting to note that the central 
(and strongly) peaked value of 7710 determined from the combined ^He and^Li 
likelihoods is at 7710 = 1.8. The corresponding value of D/H is 1.8 xlO~^, very 
close E§ to the high value of D/H in quasar absorbers 00. This is one of the 
main motivations for studying galactic chemical evolution models with high 
initial D/H and strong D destruction histories. Correspondingly, the the use 
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of reduced "^He yields is essential as described in the previous section. Since 
D and '^He are monotonic functions of rj, a prediction for 77, based on ^He and 
"Li, can be turned into a prediction for D and ^He. The corresponding 95% 
CL ranges are D/H = (4.7 - 28) x 10"^ and ^He/H = (1.3 - 2.7) x 10"^ 

If we did have full confidence in the measured value of D/H in quasar 
absorption systems, then we could perform the same statistical analysis us- 
ing ''He, "Li, and D. To include D/H, one would proceed in much the same 
way as with the other two light elements. We compute likelihood functions 
for the BBN predictions as in Eq. ( p8| ) and the likelihood function for the 
observations using D/H = (1.9 ± 0.4) x 10"'* i.e., using only the high value 
of D/H here. These are then convolved as in Eq. (^9|). In Figure 17, the 
resulting normahzed distribution, Lf^^^^{i]) is super-imposed on distributions 
appearing in Figure [T^. It is indeed startling how the three peaks, for D, ''He 
and "Li are literally on top of each other. In Figure 18, the combined distri- 
bution is shown. We now have a very clean distribution and prediction for 
V: Vw ~ 1.75j^'^ corresponding to Qsh^ = •006l'Qgj4, with the peak of the 
distribution at rjiQ — 1.75. The absence of any overlap with the high- 77 peak 
of the ^Li distribution has considerably lowered the upper limit to rj. Overall, 
the concordance limits in this case are dominated by the deuterium likelihood 
function. 

For the most part I have concentrated on the high D/H measurements in 
the likelihood analysis. If instead, we assume that the low value c3 of D/H = 
(2.5 ±0.5) X 10"^ is the primordial abundance, then we can again compare the 
likelihood distributions as in Figure 17, now substituting the low D/H value. 
As one can see from Figure 19, there is now hardly any overlap between the D 
and the ''Li distributions and essentially no overlap with the ^He distribution. 
The combined distribution shown in Figure 20 is compared with that in Figure 
18. Though one can not use this likelihood analysis to prove the correctness of 
the high D/H measurements or the incorrectness of the low D/H measurements, 
the analysis clearly shows the difference in compatibility between the two values 
of D/H and the observational determinations of ^He and "Li. To make the low 
D/H measurement compatible, one would have to argue for a shift upwards 
in ''He to a primordial value of 0.249 (a shift by 0.015) which is certainly not 
warranted at this time by the data, and a "Li depletion factor of about 3, 
which exceeds recent upper limits to the amount of depletion @. 

The implications of the resulting predictions from big bang nucleosynthesis 
on dark matter are clear. First, ii ft — 1 (as predicted by inflation), and 
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Figure 17: As in Figure 16, with the addition of the likelihood distribution for D/H. 
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Figure 18: Combined likelihood for simultaneously fitting ^He and ^Li, and D as a function 
of TJ. 
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Figure 19: As in Figure 17, with the Ukehhood distribution for low D/H. 
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Figure 20: Combined likelihood for simultaneously fitting ^He and '^Li, and low D/H as a 
function of rj. The dashed curve represents the combined distribution shown in Figure 18. 
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0.1 which is certainly a robust conclusion based on D/H, then non- 
baryonic dark matter is a necessity. Second, on the scale of small groups of 
galaxies, > 0.05, and is expected to sample the dark matter in galactic halos. 
This value is probably larger than the best estimate for ^Ib from equation (|60|). 
^Ish^ = 0.0065 corresponds to fls = 0.025 for h = 1/2. In this event, some 
non-baryonic dark matter in galactic halos is required. This conclusion is 
unchanged bv the inclusion of the high D/H measurements in QSO absorbers. 
In contrast EJ, the low D/H measurements would imply that flsh^ — 0.023 
allowing for the possibility that JIb ~ 0.1. In this case, no non-baryonic dark 
matter is required in galactic halos. However, I remind the reader that the low 
D /H is at present not consistent with either the observations of ^He nor ^Li 
and their interpretations as being primordial abundances. I will return to the 
subject of dark matter in the next lecture. 



2.^ Constraints from BBN 

Limits on particle physics beyond the standard model are mostly sensitive to 
the bounds imposed on the '*He abundance. As is well known, the *He abun- 
dance is predominantly determined by the neutron-to-proton ratio just prior 
to nucleosynthesis and is easily estimated assuming that all neutrons are incor- 
porated into "'He (see Eq. (^)). As discussed earlier, the neutron-to-proton 
ratio is fixed by its equilibrium value at the freeze-out of the weak interaction 
rates at a temperature T/ ~ 1 MeV modulo the occasional free neutron decay. 
Furthermore, freeze-out is determined by the competition between the weak 
interaction rates and the expansion rate of the Universe 

Gp'^Tf^ ^ rwk(r/) = H{Tf) ^ y/GNNTf^ (61) 

where N counts the total (equivalent) number of rclativistic particle species. 
The presence of additional neutrino flavors (or any other relativistic species) at 
the time of nucleosynthesis increases the overall energy density of the Universe 
and hence the expansion rate leading to a larger value of Tf, (n/p), and ulti- 
mateljLyp. Because of the form of Eq. ( |6l| ) it is clear that just as one can place 
limitsEJ on iV, any changes in the weak or gravitational coupling constants can 
be similarly constrained (for a recent discussion see ref. 56). In concluding 
this lecture, I will discuss the current constraint on Ni, the number of particle 
species (in neutrino units) and the limit on the strength of new interactions, if 
3 right-handed (nearly massless) neutrinos are assumed to exist. 
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In the past, ■^He (together with D) has stood out in its importance for 
BBN, because it provided a (relatively large) lower limit for the baryon-to- 
photon ratioEll, 7710 > 2.8. This limit for a long time was seen to be essential 
because it provided the only means for bounding rj from below ancLin effect 
allows one to set an upper limit on the number of neutrino flavors e3, N^, as 
well as other constraints on particle physics properties. That is, the upper 
bound to Ni, is strongly dependent on the lower bound to rj. This is easy to 
see: for lower 77, the ''He abundance drops, allowing for a larger N^, which 
would raise the ''He abundance. However, for 77 < 4 x 10^^', corresponding 
to fis/i^ < .001 — .002, which is not too different from galactic mass densities, 
there is no bound whatsoever on Of course, with the improved data on 

^Li, we do have lower bounds on rj which exceed 10^'". 

Because, of new observations of D and '^He, and the new theoretical work 
on chemical evolution sparked by these observations, the bound on N^, which 
is tied directly to these isotopes, should be called into question. As described 
earlier, the limits due to ^He are ultimately tied to the assumed yields of 
low mass stars. Using the reduced yields as depicted in Figure 14, consistent 
values oi rj < 2.8 are certainly possible. Ultimately, as I have said repeatedly, 
D/H measurements in quasar absorption systems may soon resolve this issue. 
However, the lower values of 77, relax the bounds on the number of neutrino 
flavors. 

As discussed above, the limit on N^, comes about via the change in the 
expansion rate given by the Hubble parameter, 

H-5fl,= ?^|»3M + lAiV„r (62) 

when compared to the weak interaction rates. Here iVgM refers to the standard 
model value for N. At T ~ 1 MeV, A'^sm = 43/4. Additional degrees of freedom 
will lead to an increase in the freeze-out temperature eventually leading to a 
higher ''He abundance. In fact, one can parameterize the dependence of Y on 
by 

Y = 0.2262 + 0.0131(A^ - 3) + 0.0135 Inr^m (63) 

in vicinity of ryio ~ 2. Eq. also shows the weak (log) dependence on 77. 

However, rather than use (|63| ) to obtain a limit, it is preferable to use the 
likelihood method. 

Just as ''He and ^Li were sufficient to determine a value for 77, a limit 
on can be obtained as wellll3'E3. The likelihood approach utilized above 
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can be extended to include Ny as a free parameter. Since the light clement 
abundances can be computed as functions of both 77 and N^, the likelihood 
function can be defined byE3 

LbMY. >bbn) = e-(^'-^BBN(.,iv.))V2.? (64) 



and 



L "^otai(?7, N,) = / dFiBBN {Y, Fbbn {v. N,)) Lo{Y, Yo) (65) 



Again, similar expressions are needec 
of the combined likelihood functions 



for ^Li and D. A three-dimensional view 
.21 is shown in Figure 21. In this case the 
high and low r] maxima of Figure 15, show up as peaks in the L — rj — N^, space 
(^47 when D/H is neglected and £247 when D/H is included). The peaks of the 
distribution as well as the allowed ranges of rj and N^, are more easily discerned 
in the contour plot of Figure 22 which shows the 50%, 68% and 95% confidence 
level contours in the two likelihood functions. The crosses show the location 
of the peaks of the likelihood functions. L47 peaks at Ni, = 3.0, 7710 = 1.8 (in 
agreement with our previous results E2l) and at N^, — 2.3, 7710 = 3.6. The 95% 
confidence level allows the following ranges in 77 and 

1.6 <Ny< 4.0 1.3 < 7710 < 5.0 (66) 

Note however that the ranges in r] and are strongly correlated as is evident 
in Figure 22. Since the deuterium likelihood function picks out a small range 
of values of i], largely independent of N^, its effect on ^247 is to eliminate one 
of the two peaks in L47. L247 also peaks at N^, — 3.0, 7^10 = 1-8. In this case 
the 95% contour gives the ranges 

2.0 < A^.. < 4.1 1.4 < 7710 < 2.6 (67) 

One should recall that the limit derived above is not meant for neutrinos in 
the strictest sense. That is, the limit is only useful when applied to additional 
particle degrees of freedom which necessarily do not couple to the Z°. For very 
weakly interacting particles, one must take into account the reduced abundance 
of these particles at the time of nucleosynthesiS As discussed in the first 
lecture, the number of neutrinos today is reduced relative to the number of 
photons by [T^/T^Y — 4/11. For some new particle, Xj which decoupled at 
Td> 1 MeV, the same argument based on the conservation of entropy tells us 
that 
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Figure 21: The combined two-dimensional likelihood functions for simuhancously fitting *He 
and ^Li in the top panel, and including D in the lower one as functions of both rj and A'^,^. 
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Figure 22: The top panel shows contours in the combined Ukelihood function for ^He and 
''Li. The contours represent 50% (innermost), 68% and 95% (outermost) confidence levels. 
The crosses mark the points of maximum likelihood. Also shown is the equivalent result 
when D is included. 
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Thus we can translate the bound on N^^ which is really a bound on the addi- 
tional energy density at nucleosynthesis 



^ 30 



(69) 



for additional boson states with qb degrees of freedom and fermion states with 
gp degrees of freedom. At nucleosynthesis T = = T-y and the limit Ni, < 4 
becomes 

Such a limit would allow a single additional scalar degree of freedom (which 
counts as |) such as a majoron. On the other hand, in models with right- 
handed interactions, and three right-handed neutrinos, the constraint is se- 
vere. The right-handed states must have decoupled early enough to ensure 
[Tuj^/Ty^Y < 1/3- The temperature of a decoupled state is easily determined 
from (|68| ). Three right-handed neutrinos would require N{Td) > 25, which 
from Figure 1 implies that > 140 MeV, conservatively assuming a QCD 
transition temperature of 150 MeV. If right-handed interactions are mediated 
by additional gauge interactions, associated with some scale M^/, and if we 
assume that the right handed interactions scale as M^, with a standard model- 
like coupling, then the decoupling temperature of the right handed interactions 
is related to Mz' by 

i^r-i^r (71) 

which for TdL ~ 3 MeV ( a more accurate value that the 1 MeV estimate) and 
TdL ft 140 MeV, wc find that the associated mass scale becomes AIz' ft l-G 
TeV! In general this constraint is very sensitive to the BBN limit on N^. In 
Figure 23, the allowed number of neutrino degrees of freedom are shown as a 
function of their decoupling temperature for the case oi N,y < A and Ni, < 3.3, 
shown for comparison. 



2. 5 Summary 

To summarize on the subject of big bang nucleosynthesis, I would assert that 
one can conclude that the present data on the abundances of the light element 
isotopes are consistent with the standard model of big bang nucleosynthesis. 
Using the the isotopes with the best data, ^He and ''Li, it is possible to con- 
strain the theory and obtain a best value for the baryon-to-photon ratio of 
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Figure 23: Limits on neutrino-like degrees of freedom. 

r/10 = 1-8, a corresponding value ilsh^ — .0065 and 

1.4 < 7710 < 4.2 95%CL 

.005 < ngh^ < .015 95%CL ^ ' 

For 0.4 < h < 1, we have a range .005 < fis < .10. This is a rather low value 
for the baryon density and would suggest that much of the galactic dark matter 
is non-baryonic 0. These predictions are in addition consistent with recent 
observations of D/H in quasar absorption systems which show a high value. 
Difficulty remains however, in matching the solar '^He abundance, suggesting a 
problem with our current understanding of galactic chemical evolution or the 
stellar evolution of low mass stars as they pertain to ■^He. 

3 Dark Matter 

There is increasing evidence that relative to the visible matter in the Universe, 
which is in the form of barvons. there is considerably more matter in the 
Universe that we don't seeo. Here, I will review some of the motivations for 
the existence of dark matter in the Universe. The best observational evidence 
is found on the scale of galactic halos and comes from the observed flat rotation 
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curves of galaxies. There is also mounting evidence for dark matter in elliptical 
galaxies, as well as clusters of galaxies coming from X-ray observations of 
these objects. Also, direct evidence has been obtained through the study of 
gravitational lenses. In theory, we believe there is much more matter because 
1) inflation tells us so (and there is at present no good alternative to inflation) 
and 2) our current understanding of galaxy formation only makes sense if there 
is more matter than we see. One can also make a strong case for the existence of 
non-baryonic dark matter in particular. The recurrent problem with baryonic 
dark matter is that not only is it very difficult to hide baryons, but the standard 
model of primordial nucleosynthesis would have to be discarded if all of the 
dark matter is baryonic. Fortunately, as will be covered at length in these 
proceedings, there are several attractive alternatives to baryonic dark matter. 

There is, in fact, good evidence for dark matter on the scale of galaxies 
(and their halos). Assuming that galaxies are in virial equilibrium, one expects 
that by Newton's Laws one can relate the mass at a given distance r, from the 
center of a galaxy to its rotational velocity 

M{r) oc wV/Gat (73) 

The rotational velocity, v, is measuredSH'S by observing 21 cm emission lines 
in HI regions (neutral hydrogen) beyond the point where most of the light in the 
galaxy ceases. A compilation of nearly 1000 rotation curves of spiral galaxies 
has been plotted in ref. 66 as a function of r for varying brightnesses. If the 
bulk of the mass is associated with light, then beyond the point where most 
of the light stops M would be constant and oc 1/r. This is not the case, 
as the rotation curves appear to be flat, i.e., v ^ constant outside the core of 
the galaxy. This implies that M oc r beyond the point where the light stops. 
This is one of the strongest pieces of evidence for the existence of dark matter. 
Velocity measurements indicate dark matter in elliptical galaxies as wellO. 

Galactic rotation curves are not the only observational indication for the 
existence of dark matter. X-ray emitting hot gas in elliptical galaxies also 
provides an important piece of evidence for dark matter. As an example, 
consider the large elliptical M87. The detailed profiles of the temperature and 
density of the hot X-ray emitting gas have been mapped outE3. By assuming 
hydrostatic equilibrium, these measurements allow one to determine the overall 
mass distribution in the galaxy necessary to bind the hot gas. Based on an 
isothermal model with temperature kT = 3keV (which leads to a conservative 
estimate of the total mass), Fabricant and GorensteinE3 predicted that the 
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total mass out to a radial distance of 392 Kpc is 5.7 x 10^^ Mq, whereas the 
mass in the hot gas is only 2.8 x IO^^AIq or only 5% of the total. The visible 
mass is expected to contribute only 1% of the total. The inferred value of 
based on M87 would be 0.2. 

M87 is not the only example of an elliptical galaxy in which X-ray emitting 
hot gas is observed to indicate the presence of dark matter. In addition, similar 
inferences regarding the existence of daik_matter can be made from the X-ray 
emission from small groups of galaxies eSIIJ. 

On very large scales, it is possible to get an estimate of O from the distribu- 
tion of peculiar velocities. On scale s A . where perturbations, 6, are still small, 
peculiar velocities can be expressedO as v ^ HX6Q^-^. On these scales, older 
measurements of the peculiar velocity field from the IRAS galaxv catalogue 
indicate that indeed Q is close to unity 0. Some of the new datalli indicates 
perhaps a lower value in range 0.2 - 0.5, but does not conclusively exclude 
— 1. Another piece of evidence on large scales is available from gravitational 
lensingo. The systematic lensing of the roughly 150,000 galaxies per deg^ at 
redshifts between z = 1 — 3 into arcs and arclets allow one to trace the matter 
distribution in a foreground cluster. Unfortunately, none of these observations 
reveal the identity of the dark matter. 

Theoretically, there is no lack of support for the dark matter hypothesis. 
The standard big bang model including inflation almost requires that f2 = 1 
a. If this is the case and H, — 1, then we know two things: Dark matter exists, 
since we don't see = 1 in luminous objects, and most (at least 90%) of the 
dark matter is not baryonic. The latter conclusion is a result of our previous 
discussion on BBN which constrains the baryon-to-photon ratio and hence Q b 
as in eq. ([72[). Thus 1 — fis is not only dark but also non-baryonic. 

Another important piece of theoretical evidence for dark matter comes 
from the simple fact that we are here living in a galaxy. The type of pertur- 
bations produced by inflation are, in most models, adiabatic perturbations 
(bpjp oc 5T/T)^ and I will restrict my attention to these. Indeed, the per- 
turbations produced by inflation alsoJmve the very nearly scale-free spectrum 
described by Harrison and ZeldovichLS. When produced, scale-free perturba- 
tions fall off as ^ cx I ^ (increase as the square of the wave number). At 
early times 5p/ p grows as t until the time when the horizon scale (which is 
proportional to the age of the Universe) is comparable to I. At later times, 
the growth halts (the mass contained within the volume has become smaller 
than the Jean's mass) and ^ — 5 (roughly) independent of the scale I. When 
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the Universe becomes matter dominated, the Jean's mass drops dramatically 
and growth continues as ^ cx i?- 1/T. 

The transition to matter dominance is determined by setting the energy 
densities in radiation (photons and any massless neutrinos) equal to the energy 
density in matter (baryons and any dark matter) . For three massless neutrinos 
and baryons (no dark matter), matter dominance begins at 

= 0.22mBV (74) 

and for rj <3.8x 10~^°, this corresponds to < 0.08 eV. 

Because we are considering adiabatic perturbations, there will be anisotropics 
produced in the microwave background radiation on the order of 5T/T ~ S. 
The value of 6, the amplitude of the density fluctuations at horizon crossing, 
has now been determined by COBEO, 6 = (5.7 ± 0.4) x 10"^. Without the 
existence of dark matter, dp/p in baryons could then achieve a maximum value 
of only Sp/p - AxS{T^/To) < 2 x lO"^^^^^ .^^j^ere To = 2.35 x 10""^ eV is the 
present temperature of the microwave background and A\ ~ 1 — 10 is a scale 
dependent growth factor. The overall growth in 5p/ p is too small to argue that 
growth has entered a nonlinear regime needed to explain the large value (10^) 
of Sp/p in galaxies. 

Dark matter easily remedies this dilemma in the following way. The tran- 
sition to matter dominance is determined by setting equal to each other the 
energy densities in radiation (photons and any massless neutrinos) and matter 
(baryons and any dark matter). While if we suppose that there exists dark 
matter with an abundance = n-^/n^ (the ratio of the number density of x's 
to photons) then 

T„ = 0.22tox^x (75) 

Since we can write m^Y^/GeY = D,^h'^/{A x lO"^), we have Tm/To = 2.4 x 
W^il^h^ which is to be compared with 350 in the case of baryons alone. The 
baryons, although still bound to the radiation until decoupling, now see deep 
potential wells formed by the dark matter perturbations to fall into and are 
no longer required to grow at the rate Sp/p cx R. 

With regard to dark matter and galaxy formation, all forms of dark matter 
are not equal. They can be distinguished by their relative temperature at T™ 
0. Particles which are still largely relativistic M Tm (like neutrinos or other 
particles with < 100 eV) have the property Ej that (due to free streaming) 
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they erase perturbations out to very large scales given by the Jean's mass 

Mj^Sx 10i«^^ (76) 

Thus, very large scale structures form first and galaxies are expected to frag- 
ment out later. Particles with this property are termed hot dark matter par- 
ticles. Cold particles {m-^ > 1 MeV) have the opposite behavior. Small scale 
structure forms first aggregating to form larger structures later. Neither can- 
didate is completely satisfactory when the resulting structure is compared to 
the observations. For more details, I refer the reader to reviews in refs. 62. 

Accepting the dark matter hypothesis, the first choice for a candidate 
should be something we know to exist, baryons. Though baryonic dark matter 
can not be the whole story if f2 = 1, the identity of the dark matter in galactic 
halos, which appear to contribute at the level of ~ 0.05, remains an im- 
portant question needing to be resolved. A baryon density of this magnitude 
is not excluded by nucleosynthesis. Indeed we know some of the baryons are 
dark since fl < 0.01 in the disk of the galaxy. 

It is quite difficult, however, to hide large amounts of baryonic matterEJ. 
Sites for halo baryons that have been discussed include snowballs, which tend 
to sublimate, cold hydrogen gas, which must be supported against collapse, 
and hot gas, which can be excluded by the X-ray background. Stellar objects 
(collectively called MACHOs for macroscopic compact halo objects) must ei- 
ther be so small ( M < 0.08 Mq) so as not to have begun nuclear burning 
or so massive so as to have undergone total gravitational collapse without the 
ejection of heavy elements. Intermediate mass stars are generally quite prob- 
lematic because either they are expected to still reside on the main-sequence 
today and hence would be visible, or they would have produced an excess of 
heavy elements. 

On the other hand, MACHOs are a candidate which are testable by the 
gravitational microlensing of stars in a neighboring galaxy such as the LMC 
Eil. By observing millions of stars and examining their intensity as a function 
of time, it may be possible to determine the presence of dark objects in our 
halo. It is expected that during a lensing event, a star in the LMC will have its 
intensity rise in an achromatic fashion over a periocLdi ~ 3 a/ M / .OOIMq days. 
Indeed, microlensing candidates have been found e3. For low mass objects, 
those with M < O.IM0, it appears however that the halo fraction of MACHOs 
is small, ~ 0.19j^'i[Q§i. Recent results from the MACHO collaboration show 
events with longer duration comprising about half the halo, leading to the 
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speculation of a white dwarf population in the halo. Though it is too early to 
determine the implications of these observations, they are very encouraging in 
that perhaps this issue can and will be decided. 

If we now take as given that non-baryonic dark matter is required, we are 
faced with the problem of its identity. Light neutrinos (m <_30ey) are a long- 
time standard when it comes to non-baryonic dark matter l3. Light neutrinos 
produce structure on large scales, and the natural (minimal) scale for structure 
clustering is given in Fin. ( [76[) . Hence neutrinos offer the natural possibility for 
large scale structures including filaments and voids. It seemed, however, 
that neutrinos were ruled out because they tend to produce too much large 
scale structure §1. Because the smallest non-linear structures have mass scale 
Mj and the typical galactic mass scale is ~ IO^^Mq, galaxies must fragment 



out of the larger pj 
galaxies form late 



,ce-like objects. The problem is that in such a scenario, 
^ {z < 1) whereas quasars and galaxies are seen out 
to redshifts z > 4. Recently, neutrinos are seeing somewhat of a revival in 
popularity in mixed dark matter models. 

In the standard model, the absence of a right-handed neutrino state pre- 
cludes the existence of a neutrino mass. By adding a right-handed state v^, 
it is possible to generate a Dirac mass for the neutrino, my — as for 

the charged lepton masses, where is the neutrino Yukawa coupling constant, 
and V is the Higgs expectation value. It is also possible to generate a Majorana 
mass for the neutrino when in addition to the Dirac mass term, m^iT'iiyL, a 
term MvnUn is included. In what is known as the see-saw mechanism, the two 
mass eigenstates are given by m^-^ ~ m^/M which is very light, and m^^ ^ ^'^ 
which is heavy. The state i^i is our hot dark matter candidate as V2 is in 
general not stable. 

The cosmological constraint on the mass of a light neutrino is derived from 
the overall mass density of the Universe. In general, the mass density of a light 
particle x can be expressed as 

= m^Y^n^ < Pc^ 1.06 x 10"^/io^GeV/cm^ (77) 

where Y^ ~ n^/n^ is the density of x's relative to the density of photons, for 
nho^ < 1. For neutrinos Y^, = 3/11, and one findsS 

Y^{^)my < 93eY{nho^) (78) 

where the sum runs over neutrino flavors. All particles with abundances Y 
similar to neutrinos will have a mass limit given in Eq. ([78|). 
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It was possible that neutrinos (though not any of the known flavors) could 
have had large masses, >! MeV. In that case their abundance Y is con- 
trolled by 1^,9 annihilations cil, for example, vi? ff via Z exchange. When 
the annihilations freeze-out (the annihilation rate becomes slower than the ex- 
pansion rate of the Universe), Y becomes fixed. Roughly, Y ~ (mcr^)^^ and 
p ^ o where a a is the annihilation cross-section. For neutrinos, we ex- 
nect rr 4 ~ mj^/mz^ SO that ~ \/mJ^ and we can derive a lower bound 
Eilo'E^ on the neutrino mass, > 3 — 7 GeV, depending on whether it is a 
Dirac or Majorana neutrino. Indeed, any particle with roughly a weak scale 
cross-sections will tend to give an interesting value of U,h? ^ 1. 

Due primarily to the limits from LEPo, the heavy massive neutrino has 
become simply an example and is no longer a dark matter candidate. LEP 
excludes neutrinos (with standard weak interactions) .with masses mn < 45 
GeV. Lab constraints for Dirac neutrinos are available l3, excluding neutrinos 
with masses between 10 GeV and 4.7 TeV. This is significant, since it precludes 
thii possibility of neutrino dark matter based on an asymmetry between v and 
vt3. Majorana neutrinos are excluded as dark matter since Vt^ho < 0.001 for 
nil, > 45 GeV and are thus cosmologically uninteresting. 

If we return to our example of right-handed neutrinos, we saw in the pre- 
vious section on constraints from BBN that because right-handed interactions 
are weaker than standard left-handed interactions, they decouple early and 
today are at a reduced temperature relative to i^l (cf. eq. (^)). As such, 
for Tdji > 1 MeV, n^j^/n^^ = {T,^r/T,^l)^ < 1- Thus the abundance of 
right-handed neutrinos can be written as 

Y.R^^^d){^r«^ (79) 

n.y 11 J 11 



In this case, the previous bound (78) on neutrino masses is weakened. For 
a suitably large scale for the right-haiided interactions, right-handed neutrino 
masses may be as large as a few keVtZl. 

Supersymmetric theories introduce several possible candidates. If R-parity, 
which distinguishes between "normal" matter and the supersymmetric partners 
and can be defined in terms of baryon, lepton and spin as i? = (^_i''^3B+l+2S ^ jg 
unbroken, there is at least one supersymmetric particle which must be stable. 
I will assume R-parity conservation, which is common in the minimal super- 
symmetric standard model (MSSM). R-parity is generally assumed in order to 
justify the absence of superpotential terms which can be responsible for rampid 
proton decay. The stable particle (usually called the LSP) is most probably 
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some linear combination of the i? = — 1 neutral fermions, the neutralinosEa: 
the wino W^, the partner of the 3rd component of the SU{2)l gauge boson; 
the bino, the partner of the t/(l)y gauge boson; and the two neutral Hig- 
gsinos, Hi and H2- Gluinos are expected to be heavier -rUg = (^f ) sin^ 9wM2, 
where M2 is the supersymmetry breaking SU(|2) gaugino mass- and they do 
not mix with the other states. The sneutrinoEj is also a possibility but ha£ 
been excluded as a dark matter candidate by direct t3 searches, indirect E2^ 
and accelerator searches. 

There are relatively few parameters in the minimal model and the identity 
of the LSP is determined by three: 

1 . M2 ; if one assumes a common soft supersymmetry breaking gaugino mass 
at the unification scale in the Lagrangian, C 9 M2WW . 

2. /i; a Higgs mixing mass, £ 9 fiHiH2 necessary to avoid phcnomenologi- 
cally unacceptable axions. 

3. tan/3; the ratio of Higgs expectation values, tan/3 = where vi = 
(Hi), V2 = {H2) is also a free parameter but can be chosen to be positive 
without loss of generality. In this notation, it is Hi which is responsible 
for up quark masses so that it will be natural to assume tan/3 > 1. 

The LSP can be expressed as a linear combination 

X = aW^ +1313 + 7i?i + 5H2 (80) 



LOl 



Pure state LSP possibilities are : The photinct 

^ = W^sm.ew + BcoaOw (81) 



the Higgsino, S^m 
the binJisi, 

and the Higgsinos Hy^^^ 



5° = 7?icos/3 + iJ2sin/3 (82) 
B = B (83) 



Hi2 = ^{Hi±H2) (84) 
V2 
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The solution for the coefficients a, /3, 7 and S for neutralinos that make up the 
LSP can be found by diagonaUzing the mass matrix 
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where Mi is a soft supersymmetry breaking term giving mass to the U(l) 
gaugino. In a unified theory Mi = M2 at the unification scale translates to 



5 Qi , ^ 
Ml = - — M2 

at low energies. By performing a change of basis to W^, B, A", 3° with 

i° = ^i°sin/3-^°cos/3 
the mass matrix simplifies and becomes 
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V2 
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-g2V 
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V2 -VI ,, 







f2^-f 1 



Jl — M 



B 
AO 
\ SO / 



where v'^ 



Vl 



and can be solved analj^tically. The lightest mass eigen- 



statc is the LSP. 

There are some limiting cases in which the LSP is nearly a pure state. 
When /i ^ 0, 5° is the LSP with 



2viV2 

rrig — 2— M 

When M2 ^ 0, the photino is the LSP with 

2 
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-Mo 



3(51^+522)' 

When M2 is large and M2 <C then the bino B is the LSP and 

TOS ~ Ml 



(89) 



(90) 



(91) 
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and finally when fi is large and /i <C M2 either the Higgsino state 



(12: 



or the state 



M > 



(92) 



(93) 



is the LSP depending on the sign of /i 



In Figure 24 



3, regions in the M2 , fi plane with tan (3 — 2 are shown 
in which the LSP is one of several nearly pure states, the photino, 7, the 
bino, B, a symmetric combination of the Higgsinos, ^^(12): or the Higgsino, 
S. The dashed lines show the LSP mass contours. The cross hatched regions 
correspond to parameters giving a chargino {W^, H^) state with mass rrij^ < 



AbGeV and as sue 
by LEP1.5 
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1 are exclud ed bv 



- - EI&2j. This constraint has been extended 
and LEP2 and is shown by the light shaded region 

and corresponds to regions where the chargino mass is < 80 GeV. The newer 
limit does not extend deep into the Higgsino region because of the degeneracy 
between the chargino and neutralin o. T he dark shaded region corresponds to 
an older limit on M2 from the limiil2a on the gluino mass rrig < 70 GeV or 
M2 < 22 GeV. Notice that the parameter space is dominated by the B or H12 
pure states and that the photino (often quoted as the LSP) only occupies a 
small fraction of the parameter space, as does the Higgsino combination 5°. 

The relic abundance of LSP's is determined by solving the BoltzmaniL 
equation for the LSP number density in an expanding Universe. The technique^ 
used is similar to that for computing the relic abundance of massive neutrinoO. 
The relic density depends on additional parameters in the MSSM. These in- 



clude the sfermion masses, 



the Higgs pseudo-scalar mass, uia, and the 



tri- linear masses A a s well as two phases 0^ and 6a- For binos, as was the 
case for photinos 11221, it is possible to adjust the sfermio n rn asses to obtain 



closur e de nsity. Adjusting the sfermion mixing paramet£rsll23 or CP violating 
phaseslHj'EiJ allows even greater freedom. In Figure 25Ei3, the relic abundance 
{ilh'^) is shown in the M2 — /i plane with tan (3 = 2, m j = 200 GeV. Clearly the 
MSSM offers sufficient room to solve th e d a,rk matter problem. Similar results 



have been found by other groups 



11 
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In Figur e 25. in the higgsino sector 
between -ff(i2) and 



11-; 



H12 marked off by the dashed line, co-annihilations E 
the next lightest neutralino (also a Higgsino) were not included. These tend 
to lower significantly the relic abundance in much of this sector. 
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M2(GeV) 

Figure 24: The M2-H plane in the MSSM for tan 13 = 2. 
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M2 (GeV) 



Figure 25: Relic density contours (f2/i^) in the M2 - jx plane. 
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As should be clear from figures 24 and 25, binos are a good and likely 
choice for dark matter in the MSSM. For fixed rup Q. > 0.05, for all m^j — 
20 — 250 GeV largely independent of tan (3 and the sign of /z. In addition, the 
requirement that m . > m.j^ translates into an upper bound of about 250 GeV 



1 

on the bino massli23ili3. By further adjusting the tulinear A and accounting 



for sfermion mixing this upper bound can be relaxedll23. For fixed Vlh? = 1/4, 
we would require sfermion masses of order 120 - 250 GeV for binos with masses 
in the range 20 ~ 250 GeV. The Higgsino relic density on the other hand is 
largely independent of m^. For large /x, annihilations into W and Z pairs 
dominate, while for lower /.t, it is the annihilations into Higgs scalars which 
dominate. Aside from a narrow region with ^11^12 ^ very massive 

Higgsinos with rn^j^^ < 250 GeV, the relic density of H12 is very low. Above 
about 1 TeV, these Higgsinos are also excluded. 

One can make a further reduction in the number of parameters by setting 
all of the soft scalar masses equal at the GUT scale. In particular, setting 
the Higgs soft masses equal to a common sfermion mass TOq and enforcing the 
condition of radiative electroweak symmetry breaking fixes fi and tua till- In 
this case, which we can refer to as the constrained MSSM (CMSSM) the LSP 
is almost always a binoEi3. For a given value of tan/3, the parameter space is 
best described in terms of the common gaugino masses ?ni/2 and scalar masses 
Too at the GUT scale. In Figure 26 (taken from ref. Ill) this parameter space 
is shown for tan/3 — 2. The contours labeled by a value of p represent bino 
purity contours. The other contours are labeled by the value of flh^. The 
shaded region corresponds to regions where either the LSP is not a neutralino 
(but rather a stau) or there is a light chargino. 

The TOx/2 ^ "T-o parameter space is further constrained by the recent runs 
at LEP. In Figure 27, the constraints imposed by the unsuccessful ALEPH 
chargino and neutralino scarchcslli§ (long-dashed line) at LEP 1.5 are shown. 
This search could not exclude regions of parameter space for small TOq, when 
the sneutrino mass drops bel ow the chargino mass. However, some of this area 



can be further ^^icluded c£2l by the limits on rrii} (short-dashed line), by the 
LEP limitsJlH3'E£2l on slepton production (solid line), by single-photon mea- 
surements (grey line) , and by the DO limit on the gluino mass EH3 (dot- 
ted line). The region of the plane in which 0.1 < fl^h^ < 0.3 for some 
experimentally-allowed value of /j, < is light-shaded, and the region of the 
plane in which 0.1 < ^l^h^ < 0.3 for /i determined by the CMSSM constraint 
on the scalar masses is shown dark-shaded. The constraint derived from the 
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Figure 26: The mi/2-"io plane in the CMSSM for tan/3 = 2. 
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ALEPH searches when imposing CMSSM constraints (labeled EWSB) is also 
shown as a solid line. The improvement in these bounds E2ZI by LEP 1.7 is 
illustrated in Figure 28. 

The final subject that I will cover in this lecture is the guest in n of dark 



matter detection and this topic will be covered in detail by J. Ellis The 
best chance for the detection of baryonjiiilaijk matter, as was discussed earlier, 
is through gravitational microlensingEllE^O. In principle, the duration of the 
signal should be able to determine the mass of the macho and perhaps lead to 
a determination of the composition of the galactic halo. The first test for the 
presence of neutrino dark matter is the experimental verification of neutrino 
masses. This will most likely be done first by observing neutrino oscillations 
for the case of the light neutrinos (hot dark matter). The heavy neutrinos, as 
explained earlier, are already excluded as dark matter experimentally. 

Other non baryonic dark matter detection can be separated into two basic 
methods, direct and indirect. Direct detection relies on the ability to detect 
the elastic scattering of a dark matter candidate off a nucleus in a detector. 
The detection rate will depend on the density of dark matter in the solar neigh- 
borhood, p ^ 0.3 GeV/cm^, the velocity, v ^ 300 km/s, and the elastic cross 
section, a. Spin-independent interactions are the most promising for detection. 
Dirac neutrinos have spin-independent interactions, but as noted above, these 
ve already been excluded as dark matter by direct detection experiments 
In the MSSM, it is possible for the LSP to also have spin independent in- 
teractions which are mediated by Higgs exchange. These scatterings are only 
important when the LSP is a mixed (gaugino/Higgsino) state as in the central 
regions of Figures 24 and 25. Generally, these regions have low values of Vlh^ 
(since the annihilation cross sections are also enhanced) and the parameter 
space in which the elastic cross section and relic density are large is rather 
limited. Furthermore, a significant detection rate in this case relies on a low 



_24 12^ 



mass for the Higgs scalar 

More typical of the SUSY parameter space is a LSP with spin dependent 
interactions. Elastic scatterings are primarily spin dependent whenever the 
LSP is mostly either gaugino or Higgsino. For Higgsino dark matter, Higgsinos 
with scatterings mediated by Z° avoid the -ff(i2) regions of Figures 24 and 
25, and as such are now largely excluded (the 5° region does grow at low 



tan/3E3'E23). Higgsino scatterings mediated by sfermion exchange depend on 
couplings proportional to the light quark masses and will have cross sections 
which are suppressed by [nip/mw)^, where nip is the proton mass. These 
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Figure 27: The plane with constraints from Lep 1.5. 
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Figure 28: Same as Figure 27, but updated by the LEP 1.7 run. 
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rates are generally very lowE£3. Binos, on the other hand, will have elastic 
cross sections which go as m'^/mj'^, where m is the reduced mass of the bino 
and nucleu s. T he se r ates are typically higher (reaching up to almost 0.1 events 



per kg-dayE£3'tiJ^E£il) 



Indirect methods also offer the possibility for the detection of dark matter. 
Three methods for indirect detection are often discussed. 1) 7-raysfrom dark 



matter annihilations in the galactic halo are a possible signature E£3. In the 
case of the MSSM, unless the mass of the LSP is larger than mw, the rates 
are probably too small to be detectable over background. 2) Dark matter will 
be trapped gradually in the sun, and annihilations within the sun will produce 



high energy neutrinos which may be detectedtH3; similaiiy, annihilations within 
the earth may provide a detectable neutrino signalled. This method holds 
considerable promise, as there will be a number of very large neutrino detectors 
coming on line in the future. Finally, 3) there is the possibility that halo 
annihilations into positrons and ant ipro tons in sufficient numbers to distinguish 
them from cosmic-ray backgroundsE^j. 
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